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Abstract 

The known ground state of ultrathin smectic films confined to the surface of a sphere is described 
by four +1/2 defects assembled on a great circle and a director which follows geodesic lines. 
Using a simple perturbative approach we show that for thick smectic films on a sphere with planar 
anchoring this solution breaks down, distorting the smectic layers. The instability happens when 
the bend elastic constant exceeds the anchoring strength times the radius of the inner sphere. 
Above this threshold, the formation of a periodic chevron-like structure, observed experimentally 
as well, relieves geometric frustration. We quantify the effect of thickness and curvature of 
smectic shells and provide insight into the wavelength of the observed texture. 

Keywords: geometric frustration, smectic liquid crystals, periodic structure; 


1. Introduction 

Thin liquid crystal films, coating spherical surfaces, exhibit chevron-like periodic instabil¬ 
ity patterns, as in Fig. QIIlllIl, when observed under an optical polarizing microscope. The 
mechanism of this pattern formation is two-fold. On the one hand, cooling the nematic 8CB 
(4-n-octyl-4-cyanobiphenyl) phase towards the smectic phase with the transition temperature 
TN-SmA = 33.5°C results in formation of smectic layers, orthogonal to the spherical surface, 
which locally tend to maintain constant spacing. On the other hand, because of the underly¬ 
ing curvature of a sphere and finite thickness of the liquid crystal film, this constraint cannot 
be satisfied globally, leading to geometric frustration ma. In smectic-like systems geometric 
frustration can be relieved by introducing topological defects, such as focal conics, dislocations 
or curvature wall defects Elia or forming spatially periodic texture, for example, during em- 
bryogenesis in presence of intestinal tissue growth IS]. However, a quantitative description of 
energy minimizing space-filling smectic structures still remains challenging, because of the intri¬ 
cate connection between curvature of the embedded space, local energetic constraint of equally 
spaced layers, global topological constraints as well as the boundary conditions. In this paper 
we propose a simple but realistic approach to study the effects of thickness, curvature and elastic 
properties of smectic shells on the characteristic wavelength of instability texture. 

Nematic shells with in-plane two-dimensional (2D) order have been extensively studied be¬ 
fore both experimentally and theoretically miioiin]. Smectic shells can be treated as the lim¬ 
iting case of nematic shells when the bend elastic constant is much larger than splay elastic 
constant Ki in the Frank free energy iia. The anisotropy of elastic constants + K\ influ¬ 
ences the equilibrium arrangement of topological defects in spherical nematics ifTSll as well as 
the equilibrium shape of closed vesicles with in-plane orientational order flAl . On a curved 
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Figure 1: Birefringence texture observed in experimental smectic shells (DID with radius R = 92 jim and thickness 
h = 3.4 /im. Image courtesy of Teresa Lopez-Leon and Alberto Fernandez-Nieves. 


substrate, an infinitely thin (2D) liquid crystal film with the nematic director n aligned along 
geodesics (n ■ V)n = 0 ||3|T3]|T5l contains no bend distortions, and is thus the ground state in 
the limit Ki Kj,. For a sphere the geodesics coincide with great circles, and thus the director 
n is preserved under parallel transport along meridians, resulting in splay-rich equilibrium tex¬ 
tures El. The question we pose here is how this texture changes when we allow the smectic 
film on a sphere to have a finite constant thickness. The increase of dimensionality of the system 
from 2D to 3D results in the dilation of smectic layers at the outer spherical shell, which can be 
compensated by the insertion of extra layers (dislocations) or by bending of the layers (curvature 
walls). Experimental data |[T]|2|3] suggests that geometric frustration in thick smectic shells is 
resolved by curvature walls rather than dislocations (see Fig.[^. 

In this paper we distinguish (i) ‘ultrathin’ smectic films with director field following the 
meridians of an inner sphere, which remain stable under extension into the third dimension along 
the radius and (ii) ‘thick’ smectic shells when the tilt of the director out of the tangent plane of 
the spherical surface (bending of the layers) is energetically favoured. Similar to the Barbero- 
Barberi criterion for thin nematic films ca, here we find that the configuration (i) becomes 
unstable when > RWa, where R is the radius of an inner sphere and Wa is the anchoring 
strength of the outer spherical shell. Increasing the thickness h of the smectic film above the 
critical value (oc K^jWa - R) one might expect periodic undulations of smectic layers resulting 
in the formation of crescent-like domains as in Fig.[2and in Refs. |[T]|3[3l. Thus the relationship 
between thickness h and the curvature R of smectic shell and bend elastic constant play the 
role of the magnetic field or dilative strain in the Helfrich-Hurault instability in planar geome¬ 
try UhlEl. By adopting a perturbation approach with hjR being a small parameter of our system, 
we explore the instability of the ground state (i) for h 0 with respect to the periodic solution 
in thick smectic shells. Our prediction of the critical period and the critical thickness accounting 
for chevron-like texture can be relevant to analyze experimental data and extract values of elastic 
constants and anchoring strength in the vicinity of the nematic-smectic phase transition. 

In the following we first reconsider the ground state of thin smectic shells accounting for 
splay-rich topological defects. Next, we consider the breaking down of this solution state for 
thick smectic shells, assuming infinitely strong planar anchoring at the inner shell and neglecting 
the role of splay rich topological defects {Ki <s: Kj,). In the last section we hypothesize periodic 


2 



texture of smectic layers and study the range of its stability. 

2. Ultrathin smectic shells 

The orientational order of liquid crystals is described by the director n, which is a unit vector, 
confined to tangent plane of a sphere. The in-plane spatial variations of the director held n can 
be decomposed into splay V ■ n and bend V x n elastic distortions, where V is the covariant 
derivative. Thus the elastic Frank free energy contains two independent contributions, integrated 
over surface of a sphere as 19] (TO] |TT| [l2l 

^Fei = ■ n)" + ^3|V X n|2}. (1) 

In case of nematic liquid crystals, the splay Ki and bend Kj elastic constants are of the same 
order Ki ^ Kj. Therefore both splay and bend deformations are allowed and the ground state 
of ultrathin (2D) nematic shells contains four (h- 1/2) disclinations placed at the vertices of tetra¬ 
hedron 13 . In the vicinity of nematic-smectic phase transition the elastic constant K 3 grows, 
and bend distortions are expelled from the sample, similar to the magnetic held in superconduc¬ 
tors in. The conhguration of defects in smectic liquid crystal (Kj » Ki) conhned to a sphere 
corresponds to four (h-1/2) disclinations lined up on a great circle of a sphere 113 . Below we 
illustrate this point analytically. 

Let us parametrize the director n = cos ae„ + sinaCj,, in local system of coordinates (u, v), 
with sphere of radius R being described by the radius vector r = R(cos u sin v, sin u sin v, cos v). 
The simplest director conhguration with a - njl follows meridians (or v-lines) connecting two 
H-l defects at the North (v = 0) and South (v = n) poles of a sphere as shown in Fig.[^. Thus 
smectic layers described by the level set u) - Rv - const correspond to the parallels (or u- 
lines). The director n = V«/|Vaj| = e,. is locally orthogonal to smectic layers and can be dehned 
everywhere, except at the two poles (v = 0, n), where it rotates by -i-27r. One can show explicitly 
that 


sinor a N 

V . n = —(cos V - OuCr) + — —dy,a. 


V X n = 


f?sinv 
cos a 
R sin V 


R 
sin a 

(cos V - duO) H-— o,,a, 

R 


( 2 ) 

(3) 


whence the conhguration with angle a - njl has no bend (V x n = 0). Note that, unlike the 3D 
case where the curl of a vector held is a vector (1-form), in the 2D case the curl is dehned as the 
Hodge dual of an exterior derivative applied to the 1-form, which gives a scalar (0-form) m or 
in index notation V x n = j, where is the dual of the Levi-Civita tensor ffol . Then the 
free energy 0 has purely splay contribution 


T'a=n /2 - ^ J£dudv sin vcot^ V = 2;r/ri^log^—j - ij, 


(4) 


where a is related to the size of the defect core or can be approximated by the size of 8CB 
molecule ~ 3 nm. We expect that the energy contribution of the defect core is negligible here and 
thus can be ignored. 

It was shown numerically in lfT3l and conhrmed experimentally in 0 that there is a set of 
degenerate ground states of a smectic liquid crystal on a sphere. Indeed, in the limit Kj, —» 
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Figure 2: (a) The director field n (red) on a sphere aligned along meridians (v-lines) results into two +1 disclinations at 
the North (A^) and South (5) poles. Smectic layers are orthogonal to n correspond to parallels (w-lines). (b) Splitting of 
two defects by cut-and-rotate procedure of one hemisphere, does not cost any energy da and results into four halves of 
+ 1 defect located at points A, N and B, S (other side of a sphere) separated by the angle Av = ;r/6. (c) The director field 
in (b) projected from a sphere on a complex plane with za = tan(Av/2), zb = — cot(Av/2), zn = 0 and |z 5 | = +oo. 


the bend free configuration is identical to the requirement of n being aligned along geodesic 
lines (n ■ V)n = 0 ifTSl fT4l ISll or great circles on a sphere. The geodesic on a sphere is given by 
parametric form u = M. 4 +arccos(tan cot v), depicted as blue curve in Fig.[^ By cutting a sphere 
with two +1 defects (Fig. with n = e,.) and rotating one hemisphere along a great circle by 
an arbitrary angle Av we obtain variety of states which can be transformed into one another with 
zero energy cost (see Fig. [^). The resulting director configuration on the Eastern hemisphere 
can be found by projecting two halves of +1 disclinations positioned at points za-R tan( ^) and 

Zb = -RcoK y) Fig'l};) 

0(z) = Im log [(z - Za)(z - Zb)], z^x + iy (5) 

from the complex plane z - x + iy - R tan( 5 )e'“ onto a sphere so that a(u,v) - n 12- 6(z) + arg z- 
Inserting this ansatz into 0,(0 we find Vxn = 0 and the numerical value of the integrated splay 
contribution equals to Q whence independent of Av. Thus, the ideal smectic configurations 
with equidistant layers on a sphere contain either two +1 defects at the poles (Figs. |^) or 4 
(+1/2) defects (Figs.j^) on a great circle. This result holds as long as we neglect contributions 
coming from core size of defects and describe the system as two-dimensional shell with in¬ 
plane orientational order. Below we study smectic shells with finite thickness h and analyze the 
deviation from the texture shown in Fig. |^. 

3. Finite thickness effects 

First, we extend the bend free planar solution, found in the previous section (see Fig. |^), 
into the third direction along the radius of a sphere re [R,R + h], so that we get a smectic 
shell of finite thickness h. The non-vanishing curl of the director V x e,. = e„/r now results 
in a non-uniform layer spacing in the 3D smectic shell. In the vicinity of the nematic-smectic 
phase transition both twist K 2 and bend elastic constants, entering the conventional Frank free 
energy, diverge E). For simplicity we assume that they are of the same order so /r2[ii'(Vxn)]^ + 
/r 3 [(n ■ V)n]^ = fir 3 |V X np, thus the form of the elastic free energy is analogous to Q, except for 
the region of integration. Integrating over the finite thickness gives a bend contribution InhKj, 
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and a splay contribution as 0 , multiplied by the thickness h. Indeed, the dilation of smectic 
layers extendansatzd from a spherical surface along the normal is proportional to the increase 
of the area element by dS - £/5|/,=o(l + 'IhH^ + h^K), where H - 1//? is the mean curvature 
and K - 1//?^ is the Gaussian curvature of an inner sphere. To compensate the unfavored 
increase of interlayer spacing the director n should tilt out of the tangent plane with non-zero 
gradient along the thickness, which follows from the expression for the curl of the vector field 
n = + WyC,. + HhC,, written in spherical coordinates 

Vxn = — I —[5v,(n„ sin v) - (9„n„]er H- —^—[5„nr - sin v5r(rn„)]e,,-H -[5r(rn„) - 5v,nf]e„. (6) 
r sin V r sin V r 


Let us introduce a small parameter e — h/R « 1 and the rescaled variable p 6 [0,1], such 
as r = -H ep). Then dr - ll{eR)dp ~ (9(e ') is the leading order contribution to the gradient 
of a scalar function in spherical coordinates Vcu = drOJCr + Consider a layer 

perturbation in the form a> = Rv + e^Rgip) + 0{e^), yielding 


Vcu 


V X n = 


(l +ps- + 0{e^), 

„ / 1 _ ^ 2p^-3g^-4pgr I ^ 

\R R 2R I 


(7) 

( 8 ) 


where the normal to the layers coincides with the director n, which has two non-zero components 
and accounts for the out-of-plane tilt tir oc g - dpg - g' or bend of the layers in rv-plane. We 
assume that bend deformations 0 are the dominant contribution to the free energy 0 in the 
limit » Ki (more precisely KijK^ ~ 0{e^)). The presence of pure splay topological defects 
at two poles of a sphere as well as their structure in thick smectic shells is not accounted for 
here. To find g{p) minimizing the layer dilation we solve the Euler-Lagrange equation within an 
ansatz g - p”' 


g'^ + ‘^gg'g 


+ fg" 


= 0 , 


m(6m^—7m+2)=0 


♦ gip) = Cip'^^ -I- Cap 


o2/3 


(9) 


According to experimental data IJI the inner shell p = 0 imposes strong planar alignment and 
smectic layers are not distorted .g|p=o = 0, which is also suggested by (TJ. At the outer shell 
(p = 1 ) the undulations of the director are bound by the planar anchoring, written as Ta - 
Wal2 J dS (Cr ■ n)^, which penalizes the out-of-plane tilt of the director. Up to the lowest order in 
e the total free energy is -2nR{C\ + C 2 f'E^{K^ - RWa)- The tilt g 0 is energetically favored if 
and only if > RWa- Otherwise, the in-plane alignment of the director n = (g = 0) uniform 
along the thickness is stable. The next order approximation gives the critical thickness of smectic 
shell Be oc Kil{RWa) - 1 and the amplitude of perturbation. Except for a trivial solution with 
Cl = C 2 = 0 we get the physically irrelevant cases with complex amplitudes, and two solutions 
with C 2 = 0 , given by 



Substituting this result in 0 we find up to the lowest order 

(V X n) ■ e„ = H- ^ Z 
Rrj Rt]^ 
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Figure 3: Schematic illustration of the organization of smectic layers: (a) in n’-plane described by (Tg with (counter)- 
clockwise tilt of the director corresponding to the (blue) red curves, t] = 1.05; (b) within a segment of a spherical shell. 
At the inner sphere p = 0 the director (red) is aligned along meridians, while at the outer sphere p = 1 the chevron texture 
has the highest tilt angle with non-zero projection on both u and r axes, (c) For 7] = K^KRWa) > 1 there is a critical 
wavenumber qc (extremum of the curve), characterizing an instability from the state with uniform director n and chevron 
structure with periodicity Injq. The tilt angle as well as the critical thickness grows with q. 


As predicted above and follows from ( [TT| l the tilt of the director n (|7|l, ( [TOl l, resulting in the 
bending of smectic layers in rv-plane illustrated in Fig. |^, is energetically favoured for jj > 1. 
This transition from a uniform planar state appears in the vicinity of the nematic-smectic phase 
transition accompanied by the growth of X 3 . For weak anchoring Wa at the outer shell (77 » 
1) we get V X n ~ -3epfR ~ 0(s), while for strong anchoring rj ^ 1 the tilt angle ( [TOl l is 
suppressed. Assuming the anchoring strength Wa - 10“^ J/m^ and the typical radius of a sphere 
R ^ 100 pm |[I1|2][3], one expects the critical value X 3 > 10“* J/m for the instability threshold. 
Above the threshold rj - K^jiRWa) - 1.05 and for e - hjR - 0.04 we get the estimate for the tilt 
angle as eg\p=\ ^ ±0.062. In general, the solution ( [TOl i may be homogeneous with sign, 

which is illustrated in Fig. with smectic layers bending (counter)-clockwise in polar rv-plane, 
making use explicitly of our initial ansatz 


u}^Rv±^R -^2(ep)^(^l - ^ j = const 

Note that one could have chosen or guessed the scaling of co with the power but the re¬ 
sults (|9ll-(|T^ do not change. One may guess that if the combination of both “±” and signs 
in (EH) is energetically allowed, then the change of sign of g likely results in the formation of 
curvature walls ||71. Because of the spherical geometry a naive estimate of the number of curva¬ 
ture walls is 7 r/(eg|p=i). Flowever, such a simple geometric argument based on the symmetry of 
the solution ( [T0| l may be irrelevant in accounting for the variation of birefringence colors along 
v-lines connecting two defects in Fig.E] In fact, according to CIEIEI the modulation inside cres¬ 
cent domains is interpreted as a secondary instability. In the next section we try to gain insight 
into the formation of the primary chevron-like periodic structure along the M-direction. 


ep=r/i?-l 


- = 1 ± 


|9(const-v)2 

' 8(1 - 1/7,) 


( 12 ) 


6 















4. Pattern selection mechanism 


The main goal of this paper is not to explain the complexity of fully developed texture in 
smectic shells (e.g. Fig. [T]) but to find the criterion, if there is one, for the onset of a periodic 
chevron-like solution. Again we adopt a simple perturbative approach and allow both the tilt of 
the director n in the tangent (u, v)-plane as well as variation along the thickness (r-direction). We 
assume the following form for a scalar function oj - Rv + e^Rg{p)u sin v + 0{e^), which describes 
isosurfaces of smectic layers and accounts for chevron-like texture in w-direction with period 
InRjq, where q is unknown integer (see Fig.[^). The base state with g{p) = 0 is favoured by the 
planar anchoring, but becomes unstable due to the growing bend contribution with thickness h, 
as was shown in section]^ The non-zero tilt V„(u 4^ 0 decreases the elastic energy contribution at 
the expense of the anchoring at the outer and possibly inner sphere. By optimizing the resulting 
free energy difference between the distorted state and base state, given by 


Tel + Tg 


K^R 

9q^ 


I - 4g'g''e^ + - 5g'^s^ - 12pg'g"e^ [-H 


9^2 ' 


(l+2pe)^^+6>(A (13) 


we consider the interplay of bulk elastic term and boundary term, neglecting the role of splay con¬ 
tribution. The special solution of the Euler-Lagrange equation is g(p) - ±5qp(q + p - l)/(8 Vh/r), 
linear in the number of periods of the chevron pattern q and quadratic in p, which differs 
from The integration constant Ci was determined from the natural boundary condition 

at the outer shell d{Tei + Ta)ldC\\p^\ - 0, while at the inner shell we assume g|p=o = 0. As 
before we neglected the splay /fi-term in the free energy, which may be important to estimate 
the energy of curvature walls 171. 

In Fig. the coexistence curve separates the region where the base state with unperturbed 
director n = e,. is stable from the region where the in- and out-of-plane tilt of the director with 
g{p) + 0 and ^ > 1 within assumed ansatz is energetically preferred. For q - const > 1 there 
is a critical point (qc,Sc), which identifies the wavelength InRjqc and the thickness he - ScR 
at the instability threshold. We find that for thin shells with Ec - 0.03 the chevron texture is 
fine, corresponding to large qe- For thicker shells the absolute value of the tilt angle Ecg\p=\ 
as well as the period InRjqc increases, so that the texture becomes more pronounced. In¬ 
terestingly, assuming a different perturbation ansatz in terms of the spherical harmonics a> - 
Rv + E^p"' cos{qu) sin^ v + 0{e‘^), we recover a similar prediction on the effective decrease of the 
equilibrium q with the increase of the shell thickness (e), while the exponent m 6 (j, |) depends 
stronger on q. Our results suggest the validity of the perturbation approach in the vicinity of the 
nematic-smectic phase transition and are qualitatively consistent with experimental data iniEiia 
on the formation of crescent-like domains in thicker part of smectic shells. 


5. Conclusions 

In this paper we presented a phenomenological approach to treat smectic films confined be¬ 
tween two spherical surfaces. Within the continuum Frank free energy for liquid crystals and the 
assumptions made, we have shown that the ground state of two-dimensional films with director 
aligned along geodesics of a sphere, becomes unstable for thick films when > RWa- The 
instability is driven by finite curvature of the sphere 1 jR and the requirement of the bend-free 
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configuration in the vicinity of the nematic-smectic phase transition when diverges. The re¬ 
sulting geometric frustration can be relieved by the tilt of the director out of the tangent plane of 
a sphere, as well as by in-plane undulations. The constructed simplified solution for the plausi¬ 
ble organization of smectic layers allows to identify the onset of the instability towards periodic 
chevron-like texture and provides the dependence of the critical wavelength on the critical thick¬ 
ness e - hjR of smectic shells. Future comparison between theory and experimental data □□la 
can be useful to extract the values of elastic constant Kj, and the anchoring strength Wa and even¬ 
tually propose the equilibrium 3D organizations of smectic layers within thick spherical shells. 
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